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SOLITONS IN SEMICONDUCTORS WITH A SUPERLATTICE 


A. P. Tetervov 

/1182* 

Several studies [1-4] have recently appeared which investi- 
gated the properties of semiconductors with superlattices (SL) 
in a variable electromagnetic field. The presence of a SL leads 
to very specific nonlinear effects, which are particularly apparent 
in the millimeter and submillimeter ranges. However, these 
studies devoted basic attention to calculating the nonlinear 
high frequency current. 

The propagation of strong electromagnetic waves in semi- 
conductors with SL is of definite interest. This process may be 
described by a system of Maxwell equations with a known con- 
ductivity current. The nonlinear nature of this current leads 
to the formation of nonlinear effects during the propagation. 

This article studies the possibility of nonlinear electromagnetic 
waves, namely solitons, existing in SL. 


We use the following wave equation as the initial equation 


A 


C2 



E 


4it dj 
c2 Tt ’ 


( 1 ) 


where E is the wave electric field (the field vector directed 
along the SL period parallel to the Ox axis) ; j = j (E) — the 
functional , of the field, which describes the conductivity 
current; e((o) — linear operator determining the frequency dis- 
persion of the lattice, existing in the millimeter and sub- 
millimeter ranges [5] co = co — . 


* Numbers in margin indicate pagination in foreigrt text. 
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We shall assume that the conductivity current is small as 
compared with the displacement current, i.e., the right side 
of Equation (1) contains a small parameter (the corresponding 
estimates will be given below) . Then we shall look for the 
solution of the wave equation in the form 

E = {\p^,t)exp[i((ot — kx)] + iv^,t))-^K. c. (2) 

Here (r, 0 0 exp [tcp (r, i)] • u and <}> are the amplitude and 

phase of the basic harmonics which slowly change in time and space; 
oj and k — frequency and wave vector of the waves; w{^,t) — value 

of the first order of smallness with respect to the conductivity 
current which describes the higher harmonics. 

Let us expand the current in Fourier series 

oo 

/= 2 (3) 

I~ — oo 

where Si = ! lc)t—kx+cf(rj] , and v;e place the expansion (3) 

in the wave equation. Considering the slow change in the amplitude 

and phase of the wave we obtain the following equation /1183 




w 




■ — 2tk ^\J 




I Anjj — 2ik I ijjj. 4- — Ij;,j .f _ 


4nj0) . j 

— 2 — /lexp (i(p)l exp [j (o)/ — ftx)]. 


(4) 


, o’- , a’- _ fcj 

Here ~ d? ’ — W — group velocity of the wave; 

d'^co 

f' = ^^,and the lower index of the variable \p designates 
differentiation with respect to the corresponding argument. 
Considering the slow change in the amplitude and phase of the 
wave, we may set 4'n . The component which is proportional 

to describes the dispersion of the wave responsible for the 

formation of the soliton. 
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The condition that does not contain the first harmonics 
leads to the nonlinear parabolic equation [6] , which describes the 
behavior of the complex amplitude of the basic harmonics 


i (l|;( + yri|.’:c) — Y 


Axt] + / j, (u) = 0. 


(5) 


As was shown in [3], for the quasiclassical case 
eEd, where d — the SL period; A, — width of the sth permitted 
minizone, and — does not depend on the field and energy of the 
carriers of the pulse relaxation time) in the single minizone 
approximation, the expression for jj^(E) has the form 


oo 

ii — i—\ — — ^7$t (^) Jm—l (o) j 


imcoTeft 


( 6 ) 


where /., = , - 

= ‘ edTeff 


— static current- 


r 2ne^nd-A, 


X /i 




1/2 


— generalized plasma frequency in a narrow 


zone semiconductor; n — concentration of the carriers; , — 

U / J- 

modified Bessel functions with the corresponding index; — 

Bessel function of the real argument a = p5i. At the high fre- 


quency limit, when ojTeii;:S> 


h(D 

it follows from (6) that 




B (u) 


ht 

_ Jm/t 

Lt 


a COTetf 




do (E) J I [a). 


(7) 

( 8 ) 


The quantity A (a) is the dissipative component of the variable 
current and describes the absorption of the electromagnetic wave 
by the free carriers and, as follows from the definition, 
for any type a. B(a) makes a small addition to the displacement 
current . 
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Substituting ( 7 ), (8) into the parabolic equation, we obtain 

the following system of nonlinear equations with respect to the di- 
mensionless amplitude a and the phase cfi 


(a^-) t + (a^), + 


-V Vx (a^Vi'P) 


+ iiaA (a) = 0, 


<F( + t’r9x+ Y 


+ T 


1 

2 a 


— AxQ 


Bi.a) „ 

U — — = 0 

^ 2 a 


( 9 ) 


where 


fi = 


c^k 


We assume that the solution of the nonlinear parabolic equation 
with small damping has the form of a stationary plane wave with the 
amplitude a = a^^^, and the phase 9 = 9'°’ . Let us study the 

stability of the solution obtained with respect to the small oscilla- 
tions. The solution of the linearized system ( 9 ) can be written in 
the form of waves [j(xr — 60 ]/ then 


6 = X.;; U, ± \ *4) ^ + 4tVS2^t;;x2 + ^ x2_ 


-ni/2 


( 10 ) 


where X|; and xx are the longitudinal and transverse components 


of the perturbation wave vector and Sq 



d 

(B{a)\ 

■ 1/2 

[ 2.; 

da 

1 “ I. 



From an analysis of the expression obtained, we find that at 

y'<0 and ^^<0 a plane wave is unstable with respect to the 
longitudinal perturbations with 5Cx = 0 , and small values of 
x2|]<4|55l . As is known, this instability leads to self- 
modulation of the wave (its decomposition into solitons) [6] . 

If , then from (8) 

, [i sin 2a<°> 

’ ( 11 ) 

and then a plane wave with the amplitude , 

1), m = 0,1,2,... ^^2) 
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decomposes into solitons. We should note that the condition 
corresponds to a very strong field 

ed 


In the case (or, which is the same thing, the wave 

field is very weak, i.e., ), the sign of Sq, as can be 

readily shown, is determined by the sign of v^. At y'>o '*^he 

solution in the form of a plane wave is stable and at the 

instability leads to self-modulation of the wave. 


The sign of K can be determined from the dispersion equation 
of linear theory, considering the lattice as a set of linear 
oscillators with an eigenfrequency of and with the force 


k-c^- 


= e(o)) = eo 



(13) 


where 


is the dielectric constant of the lattice, 

Q 8 (£0) 

V i— L 0 


At CO- > co2 


(14) 
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and at co^ coj 


— 2 


CO- 


o 


co?ft'^ W 


•< 0 . 


(14a) 


To determine the form of the soliton, we transformed the 
system (9) by changing \=-x — vrt, i = vj. to the form 


(a2), + (a2(p.)g+-^fl/l (a) = 0, 
r 

,1 , B(a) , ^ 


(15) 


Disregarding the damping, the solution (15) can be sought 
as a function of ^1 = 1 — i.e., in the form of stationary waves. 

Here W is the "velocity" of a stationary wave, determined from the 
boundary conditions. It follows from the upper equation of the system 
(15) 

' u- 

(16) 


a~ — a^ 

= r ' 
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where is the integration constant. Substituting the expression 
obtained in the lower equation of System (15) , we obtain the 
equation of the nonlinear oscillator: 




- 117 - 



— 4 s (o) = 0 . 


(17) 


Multiplying it by 2a„ and integrating, considering (8) we obtain 




(18) 


where a .2 is a constant. The solution of Equation (18) can be 
represented in the form 


a(n) 

,1 


,1 

1 

da 


1 

ly 



+ 

1 1 

1/2 


(19) 


(aQ is the initial value of the wave amplitude) . 

If a increases with an increase in n , we must select the 
upper sign in (19) . If it decreases, we select the lower sign. 
The study of the forms of stationary waves is similar to that 
of [7]. 


Let us consider the case of a weak field «<1. Expanding the 
Bessel function in terms of the small argument, we reduce (19) to the 
form 


<i(n) 


da 


a, 4- 




lt/2 


( 20 ) 
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At 02 = — ^, a, = 0 


(20) is readily integrated, as a 


result of which at ^^r<0 we obtain the soliton equation 


4 / o I y” I \ r 


2\v' 


1/2 n 


( 21 ) 


The more general solutions of the equation (20) describe the stationary 
periodic waves and can be expressed in terms of elliptical functions. 


There is a solution in the form of a single wave for large 
values of the wave amplitude. To prove this, we note that when the 


conditions 


a — a, 


1 and ao>l are satisfied (18) assumes the 


form of the Equation of sin-Gordon [3] 

H sin2 a = C {a^). 

1 OTQo V°' 


( 22 ) 


Here C(aQ) is a constant. The solution of Equation (22) has the 
form 

a(r\) 

da 

Tl/2 • 

(23) 




a(r\) 


^ (^o) sin^ a 


a 


In particular, assuming = C'(ap) = 0, 

we obtain the single equation 

a (ii) = 2arctg I 


1 ^ 

/ n 1 1)' 1 Qp \ 

exp 1 ± j 



(24) 


(the derivative a has the usual form of the soliton: 


= ± ( — jt: — 1 sech 




(24a) 


It can be readily shown that the soliton solutions (21) and 
(24) are stable with respect to small perturbations. 



In conclusion, we shall make certain numerical estimates. The 
condition for the smallness of the conductivity current as compared 
with the displacement current means the following 

CD=eo»-^' + (25) 

Thus, the conductivity current in the semiconductor with eo~ 10, 10‘® cm~^ 

and Ai?s;10-2 is small as compared with the displacement 
current at frequencies of co»5-10'2 sec^. In the region of 
these frequencies, the condition a = 1 corresponds to the field 
strength of the wave E = 3000 V/cm, 

The author would like to thank F.G. Bass for his stimulating 
discussions. 
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